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ABSTRACT 


The main purpose of this study is to gain insight into physics of the 
continuum spectrum of active galactic nuclei (AGNs) using a large data set and 
rigorous statistical methods. We have constructed a database for 469 objects 
which include radio selected quasars, optically selected quasars, X-ray 
selected AGNs, BL Lac objects and optically unidentified compact radio 
sources. Each object has measurements of its radio, optical, X-ray core 
continuum luminosity, though many of them are upper limits. Since many radio 
sources have extended components, we carefully select out the core component 
from the total radio luminosity. With 'survival analysis' statistical 
methods, which can treat upper limits correctly, these data can yield better 
statistical results than those previously obtained. 

A variety of statistical tests are preformed, such as the comparison of 
the luminosity functions in different subsmples, and linear regressions of 
luminosities in different bands. Interpretation of the results leads to the 
following tentative conclusions: (1) The main emission mechanism of optically 
selected quasars and X-ray selected AGNs is thermal, while that of BL Lac 
objects is synchrotron; (2) radio selected quasars may have two different 
emission mechanisms in the X-ray band; (3) BL Lac objects appear to be special 
cases of the radio selected quasars; (4) some compact radio sources show the 
possibility of SSC in the optical band; and (5) the spectral index between the 
optical and the X-ray bands depends on the optical luminosity. 


I. INTRODUCTION 


One of the most important problems in the studies of active galactic 
nuclei (AGNs) is understanding the mechanisms of underlying continuum 
emission. Although there are already very many observations of AGNs across 
the whole range of spectrum, and the knowledge of properties of AGNs has been 
improving considerably, we still do not understand the fundamental emission 
mechanisms. The continuum emission spectrum distributions probably arise from 
combinations of several mechanisms, including both thermal and non-thermal 
processes. The radio emission is thought to be incoherent synchrotron 
radiation. For some objects like BL Lac objects, the synchrotron spectrum 
clearly extends to optical region and perhaps to the X-ray band. On the other 
hand, most of optically selected quasars do not show radio emission and have 
unpolarized continua. The infrared to ultraviolet regions often have spectral 
bumps that are not well understood. The emission mechanisms may be 
combinations of unpolarized synchrotron, bremsstrahlung from an accretion 
disk, dust emission, stellar photospheric emission, or Compton scattering by 
thermal or non-thermal electrons. In the X-ray band, synchrotron self-Compton 
(SSC) is one of the more popular models, though multi -temperature 
bremsstruhlung is also probable. 

Since the Einstein Observatory has provided high quality X-ray 
observations, statistical studies of AGN continua are flourishing (Ku, 

Helfand, and Lucy 1980; Zamorani et al . 1981; Owen, Helfand, and Spngler 1981; 
Owen and Puschell 1982; Kriss and Canizares 1982; Reichert et . al. 1982; 
Zamorani 1982; Avni and Tananbaum 1982; Blumenthal, Keel, and Miller 1982; 
Kembhavi and Fabian 1982; Schwartz and Ku 1983; Ledden and O'Dell 1983; 
Tananbaum, Wandle, and Zamorani 1983; Katgert, Thuan, and Windhorst 1983; 
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Marshall et al. 1983; Bregman 1984; Maccacaro et al . 1984; Henriksen, 

Marshall, and Mushotzky 1984; Cruz -Gonzales and Huchra 1984; Miller 1984; 
Marshall et al . 1984; Ledden and O'Dell 1985; Kriss and Canizares 1985; Stocke 
et al . 1985; Franceschini , Gioia, and Maccacaro 1986). Considerable attention 
has been focused on the evaluation and interpretation of the average optical - 
to-X-ray spectral index, <a QX > > f° r various samples of AGNs. These finding 
are briefly summarized in Table 1. The < a ox > values have been used, for 
example, to infer that radio -selected quasars are several times more X-ray 
luminous than optically selected quasars (Ku, Helfand, and Lucy 1980, Zamorani 
et al. 1981), that the broad band spectral index evolves with redshift in 
radio-quiet quasars (Ku, Helfand, and Lucy 1980, Zamorani et al . 1981), and 
that high polarization quasars and BL Lac objects have similar continuum 
shapes (Ledden and O'Dell 1985). 

In addition to the comparison of broad and spectral indices, statistical 
correlations between radio, optical, and X-ray emissions in AGNs have also 
been studied by the previous workers. Ku, Helfand, and Lucy (1980) show a 
general correlation between radio and X-ray emissions in quasars, which was 
confirmed and refined in our examination of radio-loud quasars (Kembhavi, 
Feigelson, and Singh 1986; hereafter Paper II). The close correlation between 
optical and X-ray luminosities in Seyferts and quasars has been established by 
a number of workers (e.g. Reichert et al. 1982, Blumenthal, Keel, and Miller 
1982, Kriss and Canizares 1983, Kriss and Canizares 1985). Zamorani (1984) 
and Avni and Tananbaum (1986) have examined the relations between the spectral 
index a and optical luminosity and redshift, and Zamorani (1984) has raised 
the possibility that the quasars X-ray luminosities are simultaneously 
correlated with their radio and optical luminosities. 
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Some researchers have expressed a reluctance to examine directly the 

correlations between luminosities at different spectral bands for fear of 

encountering spurious correlations, as all luminosities for a given object are 
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scaled by the same ( distance ) factor. We have shown, however, that no such 
distance -dependent effect occurs if all objects are considered, including 
those not detected (Feigelson and Berg 1983 and Paper II). A major reason 
this study has been undertaken is that powerful and well-established 
statistical techniques are now available that fully account for the presence 
of upper limits in luminosity- luminosity diagrams (Feigelson and Nelson 1985, 
Isobe, Feigelson, and Nelson 1986). 

The present study represents an improvement upon previous studies in 
three respects. First, following Paper II, we take particular care to 
consider CORE radio emission rather than TOTAL radio emission from each AGNs. 
The total radio flux frequently includes jets and lobes, which does not 
reflect the current state of activity in the nucleus. While this distinction 
is small for some classes of AGNs (e.g. radio selected BL Lac objects), it is 
a considerable correction for others (e.g. 3C and 4C quasars). Note that our 
earlier VLA observations (Feigelson, Isobe, and Kembhavi 1984; hereafter Paper 
I) were specifically designed to acquire radio core fluxes for this study. 
Radio observations from certain optically selected quasars not presented in 
Paper I, are now presented in § II. Second, we analyze a much larger number 
of objects than earlier studies by virtue of having collected most or all of 
the extant literature. The database (§ III) includes all AGNs (except for 
certain classes, such as Seyfert galaxies and radio galaxies, where the data 
suffer significant ambiguities) for which radio core, optical and X-ray 
observations have been reported. All upper limits are included. Third, we 
use the wide variety of statistical methods provided by 'survival analysis', 
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the field of applied statistics developed over several decades to solve 
problems involving upper limits in medical and industrial situations. These 
methods are reviewed in § IV; the reader is encouraged to examine Feigelson 
and Nelson (1985) and Isobe, Feigelson, and Nelson (1986) for more details. 
Applying these statistical methods to the database in § III, we calculate the 
correlations and linear regressions between X-ray, optical, and radio 
luminosities . We also investigate several specific issues: (1) The dependence 
of optical to X-ray spectral index on the optical luminosity and the redshift; 
(2) a proposed two - component model for X-ray emission of the radio selected 
quasars; (3) a comparison of the BL Lac objects with the radio selected 
quasars; and (4) a comparison of the optically selected quasars with the X-ray 
selected AGNs . The results from these investigation and their interpretations 
are presented in § V. § VI summarizes the whole study. 
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II. RADIO OBSERVATIONS 


Although most of the data are drawn from the published literatures, we 
have made two sets of observations with the NRAO^ Very Large Array (VLA) to 
improve the quality of radio data on certain AGNs with measured X-ray 
luminosities . 


The National Radio Astronomy Observatory is operated by Associated 
Universities, Inc., under contract with the National Science Foundation. 


On 23-24 October 1982, thirty- six optically selected quasars with X-ray 
properties measured by Ku, Helfand, and Lucy (1980), or Zamorani et al . (1981) 
were observed at the VLA. These radio quiet quasars were observed at the same 
time as the radio loud quasars discussed in Paper I. The array was in the 
standard B configuration and 26 antennas were operating. The results are 
given in Table 2a. Twenty-eight quasars were not detected with 5 x rms upper 
limits around 1 mJy, and 8 quasars were detected with flux densities between 
0.8 and 40.5 mJy. Of these, four were not (to our knowledge) previously known 
radio sources, including the comparatively bright quasars GQ Com and V396 Her. 

In the second set of VLA observations, we observed the X-ray selected 
AGNs from the serendipitous Einstien IPC sample of Kriss and Canizares (1982) . 
The observations were performed along with the survey of discussed by Gioia 
et al. (1984) on 28-30 November 1981 with the VLA in the C configuration. The 
results are given in Table 2b. Data for one object in the sample, 0514-003, 
were not good. Snapshots of -12 minutes duration gave 5 x rms upper limits 
around 0.7 mJy for 21 of the sources. Three are detected, one of which 
(1401+085, z=0.43, S^=-18.8 mJy) is quite radio luminous. 
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III. DATABASE 


This study is base on the radio, optical and X-ray luminosities of a 
variety of AGNs given in Tables 3, 4, and 5. In Table 3, we show data on the 
radio selected, optically selected, and X-ray selected samples of emission 
line AGNs. Data on BL Lacertae type objects are shown in Tables 4a and 4b. 
Table 5 includes optically faint or undetected AGNs for which redshifts 
measurements are not available. 

The database, however, excludes certain classes of AGNs for which 
unambiguous radio cores, optical magnitudes or X-ray data are not available. 
Few Seyfert galaxies have optical core magnitudes reported separately from the 
host galaxy brightness, and their radio structures are frequently complex so 
that the core is not readily discriminated from ejecta (e.g. Ulvestad and 
Wilson 1983). Radio galaxy nucleus X-ray emission may be often confused by 
the diffuse X-ray of the surrounding intercluster medium (Feigelson and Berg 
1983), and again optical core magnitudes are usually not available. The PG 
sample of bright optically selected quasars (Schmitt and Green 1986, Tananbaum 
et al . 1986), the Braccesi and other fields of faint optically selected 
quasars (Braccesi et al. 1970, Marshall et al. 1984) have X-ray observations, 
but sensitive radio measurements have yet to be published. All of these 
samples have therefore been omitted from our study. 

Data in Tables 3 and 4 are organized as follows: Column 1 lists the 
source by its Right Ascension and Declination. Column 2 gives the catalog 
name from the various radio and optical surveys. Column 3 gives the redshift 
values taken from the X-ray literature, if available, or from other sources 
described in the Notes to Tables 3 and 4. In column 4, we give the radio core 
luminosities for the sources. The luminosities are computed using the 
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following formula: 


i- 4 it f (l+z) (c * X) 


where a is the spectral index within the appropriate spectral band (f-i/ ) , 

and f is the observed flux density, z is the redshift and is the luminosity 


distance given by 


cz . 1+2 

d r ^ < 1 + t _ > • 


where we assume a Hubble constant H q- 50 Mpc/km/sec and q^O. For the radio 
core emission, we assume the spectral index within the radio band, a -0. 

Since many radio selected quasars have extended components, we use the 
following procedure to find the core luminosity density at 5 GHz: 

i) If a map that clearly resolves the core from any jets or lobes is 
available, the core flux density is used. If the map is not at 5 GHz, the 
core spectral index is assumed to be 0.0. 

ii) If the source is fully resolved and the core is not detected, we use an 
upper limit given in the literature. If an upper limit is not given, the flux 
density of the weakest component is used as an upper limit. 

iii) If an interferometric map is not available, but the source is seen by 
single dishes and the spectral index between 1.4 GHz and 5 GHz is less than 
0.3, the entire flux density given for the object is assumed to be the core 
flux density. 

iv) If single dish data are available, and the spectral index is steeper than 
0.3, the 5 GHz flux density given for the object is treated as an upper limit, 
even if it is detected. This is because of the probable existence of extended 
components . Although a distribution of the upper limits set by this procedure 
may not be same as that of the upper limits due to the flux limited 
observation, we assume that all upper limits belong to a same population. As 
discussed in Paper II, even such careful efforts to isolate radio core fluxes 


9 



can leave . a residual - 50% extra flux from VLBI scale jets. 


The majority of optically selected quasars and X-ray selected AGNs are 

not detected in the radio band. Those that are detected are generally faint 

and unresolved, and their flux is assumed to arise from a core with a =0.0. 

r 

Cases where the radio measurement, either detection or upper limit, was made 
at 2.7 GHz or 1.4 GHz rather than 5 GHz are marked by a * or + in Table 3. 


The optical luminosities are given in column 5. For optical emission, 
the spectral index within the optical band, a^-1.0, is assumed. Visual 
magnitudes are mainly taken from Hewitt and Burbidge (1980) . The visual 
magnitudes are converted to the optical luminosity density at 2500 A according 
to Zamorani et al . (1981), 


0.072 


where V 


corr 


log( V ~37.878 + 2.1og[z(l + f)] - 0.47^ + 
is the visual magnitude corrected for Mgll line emission (Schmidt 


1968), and the last term is a correction for galactic absorption, where b is 


the galactic latitude. If only a blue magnitude is available, the equation is 
modified as 


log(i o )-38.0U 4 2>log[z(l+f)] - 0.4 4 ± 221 , 

where B corr is the corrected blue magnitude from Schmidt. If a redshift is 

not available, an optical flux density at 2500 A at the observer's frame is 

computed to find a and a , 
ro ox 

log (f )=-19 . 756 - 0.4V + - °; ~r . 

° o sin(b) 

For some compact radio objects, since only red magnitudes are available, we 
need to change the constant in the last equation to -19.521 (Johnson 1966). 


Column 6 lists the X-ray luminosities computed for the 0. 5-4.5 keV energy 
range in the emitting frame and assuming the spectral index within the X-ray 
band =0.5. The X-ray data are obtained mostly from the observations with 
the Einstein Observatory , although some data from observations with the HERO-1 
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satellite are also included. If the X-ray luminosity in the above energy band 
already exists, then it is adopted directly, else it is calculated from the 
observed flux in the 0.5-4. 5 keV energy band according to the fomulae given 
after Table 5 . Column 7 lists the radio spectral index of the core wherever 
the measurement exists. 

The values for the spectral index ( a ) computed between the radio (5 
GHz) and the optical (2500 A) bands in the emitting frame are listed in 
column 8. These values are calculated using the following expression: 

a -(logi - logi )/5.38. 

If the redshifts are not available, we use the flux densities instead of the 
luminosities. The spectral index (a ) between 2500 A and 2 keV emission is 
given by 

a -(logi - logi - 17 . 98)/2 . 61 . 
ox e o & x ' 

The values for a Qx are listed in column 9. Column 10, 11, and 12 list the 
references for the radio, optical, and X-ray data respectively. 


Although not shown as separate tables, a few other subsamples are used. 
For some statistical problems, we use spatially resolved radio selected 
quasars, unresolved radio selected quasars with flat (a^<0.3), and steep 
(a^>0.3) spectra. These samples are discussed in detail in Paper II. 


Comparison of our data calculated luminosities, a and a values to 

ro ox 

previous collection of continuum emission in AGNs , such as Ku, Helfand, and 
Lucy (1980), Zamorani et al. (1981), and Ledden and O'Dell (1985), shows 
relatively good agreement. One relatively large discrepancy is the radio 
luminosity, since we use the core instead of the total luminosities. The 
difference in the radio luminosities often reachs an order of magnitude 
difference. Because of this, the spectral index between the radio and optical 


band can show relatively large discrepancies . Ii 


addition, the BL Lac ouj 6 C ls 
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often show large differences from other studies . This may be due to their 
high degree of variability. 

There are several possible sources of error in the database. First, data 
are collected from variety of references and may be differently treated in 
each reference. Second, and related to the first point, uncertainty arises 
when a source is variable and observations are not done simultaneously. 

Third, some error is caused by the extrapolation of published data with fixed 
spectral indices to compute luminosity densities at consistent wavelength in 
the emitting frame. For example, this may be an important error source for 
the optical luminosity density, since we do not consider the effect of 3000 A 
bump and other effects in UV region. Fourth, although we use point radio 
sources in VLA maps as cores , VLBI observations often show that the cores can 
be resolved to further small scale. Therefore, the core radio luminosity 
density may be systematically overestimated. Some small errors also arise 
from assuming a specific cosmology. Although the Hubble constant H ^ does not 
affect correlations, the cosmological constant q ^ does. 

Based on comparisons of luminosities with previous studies and our 
estimation of the size of these possible sources of error, we find typical la 
uncertainties of ±0.2 in log(i^) , log(i Q ) , and logCi^) . Although these error 
sources may seem large, uncertainties of less than 0.5 in log form are not 
very significant, since the ranges of 2^, 2 q , and 2 are frequently 10^. 
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IV.- STATISTICAL METHODS 


Since the data set contains many upper limits, survival analysis must be 
used to treat the data correctly. Survival analysis or the analysis of 
lifetime data, is developed over several decades to deal with problems arises 
in clinical epidemiology, actuarial science and industrial reliability, where 
'censored data' (i.e. upper or lower limits) frequently arise. The methods 
are typically extensions of parametric (e.g. least square regression) or non- 
parametric (e.g. Kolmogorov- Smirnov or Mann-Whitney tests) statistical tools 
used for uncensored data, and frequently involve maximum- likelihood concepts. 
Most of the procedures we use in this study are described by Feigelson and 
Nelson (1985) and Isobe, Feigelson, and Nelson (1986). The former study 
treats problems involving one variable: the Kaplan-Meier estimator is the 
maximum- likelihood estimator of the luminosity and gives a mean luminosity and 
a standard deviation for a sample; the Gehan and logrank tests measure whether 
two subsamples are drawn from a same parent population. The latter study 
treats correlation and regression between two variables: Cox regression and 
the generalized Kendall's r test (the BHK method) which measure the degree of 
independence; the EM algorithm and Buckley- James methods perform linear 
regression on the data. 

One new method is used in this study. Previously, we could not fit a line 

on a data set which contains upper (or lower) limits in both independent and 

dependent variables, except by Schmitt's (1985) method which does not provide 

analytic estimates of the uncertainties for the regression parameters. Using 

the BHK method described by Isobe, Feigelson, and Nelson (1986), we have 

developed a method to find a slope coefficient and uncertainty. Consider a 

database in two variables (X . Y .) with possible non- detections in both 

-*• . 1 
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variables. For a range of slope coefficients b, calculate residuals r -X -bX . 
The value of b that minimizes the generalized Kendall's r rank correlation 
coefficient between the r and X is the most probable value. The la 
uncertainties can be obtained by finding the slope coefficients which give 31% 
of the maximum probability. To find the intercept coefficient, we use the 
Kaplan-Meier estimator. First, get the residuals r with the best slope 
coefficient b. The best estimate of the intercept coefficient is the Kaplan- 
Meier mean of the residuals. This combination of survival analysis methods on 
doubly censored data may permit parameter estimation for non-linear models as 
well. We use it in § IV to test a two-component model of quasar X-ray 
emission. Although similar procedures have been already suggested for non- 
censored data sets by statisticians (Sen 1968, Efron 1984, Lancaster and Quade 
1985), there is no statistical study for censored data sets yet. From our 
experience, however, the resulting regression coefficients are quite 
satisfactory when compared to those obtained by other methods. 

Using these survival analysis techniques, we analyze our database. Cox 
regression and the BHK method are used to compute the correlation 
probabilities between the radio, optical, and X-ray luminosities, and the EM 
algorithm and Buckley-James method (and the new linear regression method, if 
needed) are used to calculate the linear regression coefficients (§ Va) . The 
mean values of the spectral indices are calculated by Kaplan-Meier estimator 
(§ Vb) . Multi-dimensional linear regression among the optical luminosity, the 
redshift, and the spectral index between the optical and the X-ray bands is 
discussed in § Vc. The new regression method is applied to analyze the two- 
component model for the X-ray emission of the radio selected quasars (§ Vc) . 
For the comparison of the BL Lac objects and the radio selected quasars, and 
the comparison of the optically selected and the X-ray selected AGNs , the two 
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sample tests (Gehan and logrank tests) are used (§ Ve,f). 
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V. STATISTICAL ANALYSIS AND INTERPRETATION 


We now proceed to investigate a number of statistical relations between 

the Z , Z , and Z values listed for the various samples in Tables 3 to 5 . 
r o x r 

Most of the relationships are illustrated in Figures 1 to 3 , which plot the 
luminosity densities against each other, and Figure 4, which plots the 
interband spectral indices a vs. . The plots are displayed so that the 
various subsamples can be easily distinguished. 

All data lie in the range of 29 < log(i^) <37, 28 < log(i o ) <33, 42 < 

log (Z ) < 48 , -0.3<a <1.2, and 0.8 < a <1.9. The radio quasars tend 

x ro ox 

to occupy higher and relatively wider range (six decades for radio, four 
decades for optical and X-ray). BL Lac objects occupy similar range as the 
radio quasars. The optically selected quasars are usually one order of 
magnitude weaker than the radio quasars. The X-ray selected AGNs occupy same 
range as the optically selected quasars do. The X-ray selected BL Lac objects 
are few in number and occupy only two decades in any luminosity ; hence we will 
not be surprised if significant statistical results are not obtained. 

a) Correlations and Linear Regressions between AGN Luminosities 

Using Cox regression and the BHK method, we establish the significance 

level of correlations between radio and X-ray luminosities, optical and X-ray 

luminosities, optical and radio luminosities, and a and a for all 

ro ox 

subsamples described in § III. Quantitative results are shown in Table 6a. 

The first column lists the name of the samples and the second column lists the 
correlations tested. The third column shows the total number of the objects 
and the fourth column lists the number of censored data. Three values are the 
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numbers of data points which are censored in the independent variable only, 
the dependent variable only, and both the variables, respectively. The fifth 
and sixth columns show the correlation probabilities by Cox regression and the 
BHK method; This is the probabilities that the two variables are not 
correlated with each other. The last column identifies the corresponding 
figure. Except for some spectral index correlations, all subsamples show high 
significant level between all spectral bands. For example, even the optically 
selected quasars show a highly significant correlation (P<0.01%) for the radio 
and X-ray luminosity relation. 

Although we find very high significant levels for nearly all 
correlations, it is difficult to tell which correlations are intrinsic and 
which correlations are secondary. For a completely detected data set, we can 
use a partial linear correlation method and a partial rank correlation, but 
these methods cannot treat a censored data set. Using the generalized 
Kendall's r correlation coefficient, we may be allowed to use a normal partial 
rank correlation formula, but since a partial rank correlation is distribution 
free, we cannot get significant levels. We show the partial correlation 
coefficients results shown in Table 7 , but they can be used only for 
qualitative examintion. 

For the radio selected quasars, the optical / X-ray relation is most 
significant (r-0.45, where r is the partial correlation coefficient in Table 
7) and the radio / X-ray relation is also important (r-0.42), but the optical 
/ radio relation may not be significant (r-0.13). For the optically selected 
quasars and the X-ray selected AGNs , the optical / X-ray relation is most 
important and two other relations may not be significant. For the radio 
selected BL Lac objects, the radio / X-ray relation is most significant and 
the optical / X-ray relation is moderately significant, but the optical / 
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radio relation is weaker. For the X-ray selected BL Lac objects, the optical 

/ X-ray relation is most important, and the optical / radio relation is 

moderately important, though because of small size of the data, this finding 

may not be accurate . We thus find that the Z / Z relations are typically of 

greatest importance , with the Z / Z and Z / Z relations important only in 

r x ir o 

certain subsamples . 

The linear regressions are done mainly by the EM algorithm which assumes 
the luminosities are distributed in a Gaussian distribution about the best fit 
line, and the Bukley- James method, which makes no assumptions regarding the 
distribution of residuals. Since the relation between log(i^) and log(^) 
contains upper limits in both variables in the same subsamples, the new method 
described in § IV is used to compute coefficients. Quantitative results are 
shown in Table 6b. The first column lists the name of the samples and the 
second column lists the independent and the dependent variables. The third 
column shows the total number of the objects and the forth column lists the 
number of censored data. The fifth and sixth columns show the linear 
regression results by the EM algorithm and Buckley- James method respectively. 
The first row in the each set shows the intercept coefficient, the second row 
shows the slope coefficient, and the last row shows the standard deviation. 

If only one set of the result appears, the regression was done by the new 
method described in the § IV, or by a normal least square method, if there are 
no censored data. For example, the radio selected quasars have the linear 
regression form, log(^ x )=29 . 0+(0 . 48±0 .OSJlogCi^) . This best fit line is shown 
in Figure la. 

The linear regressions can be summarized as follows. The radio/X-ray 

correlation is about for quasars of all types but is significantly 

0 8 

steeper (i «r ' ) for BL Lac objects. The optical/X-ray correlation behaves 
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for BL Lacs) , though the 


similarly ^ for quasars compared to 

X-ray selected AGN subsample does not fit the pattern. The radio/optical 
correlation is for the radio and the optical quasars as well as the BL 

Lac objects, but, as stated above, may be an indirect consequence of the 2 2 
and % 0 /& x correlations. The correlation is present in X-ray selected AGNs and 
the X-ray selected BL Lac objects, but they have very different forms from the 
others . 

A plausible theoretical interpretation of these results might be as 

follow. First, the results of partial correlation analysis (strong 

correlation between 2 and 2 but not between 2 and 2 ) suggest that for the 

ox r x 

optically selected quasars and the X-ray selected AGNs, the thermal emission 
(e.g. bremsstrahlung, Comptonization) is the dominant mechanism. The thermal 
emission scale according to ^ . This form is predicted by Tucker 

(1983) and Schlosman, Shaham, and Shaviv (1984) for the thermal emission from 
accretion disks. Second, for the BL Lac objects, because of the strong 
correlation between 2 ^ and 2^, non- thermal emission is most important. The 
non- thermal emission may have the form of ^ • Third, for the radio 

selected quasars, the X-ray emission depends on both the radio and optical 
emissions. The result indicates that there are possibly two different X-ray 
emission mechanisms (see § Vd for the further disccusion) . 

b) Mean Values of the Spectral Indices 


The relation between the interband spectral indices are shown in Figures 
4 a to 4e with quantitative results given in Table 8. All mean values and 
standard deviations were calculated using the Kaplan-Meier estimator (§ IV). 
As expected from simple selection effects, <0c ro > is relatively large in the 
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radio selected quasars and BL Lac objects, and <<* ox > is relatively large in 

the optically selected quasars . All subsamples except the compact radio 

sources show <a > - 0.5 and <a > > 1.0. On the other hand, the compact 
ro ox > r 

radio sources show <a > - <a > - 1.0. This means that one power law with 

ro ox r 

a-1.0 can express the entire emission between the radio and X-ray bands for 
this subsample. Close examination of Table 5 also tells that some objects may 
have <a > > <a >. Since this relation cannot be readily produced by thermal 
or the synchrotron processes, this may be direct evidence for the SSC model 
operating in the optical to X-ray bands. 

In Table 1, we summarized published values a and a from recent 

ox ro 

literature . A comparison with our results shows that although they are not 
exactly the same as our values , they agree reasonably well . For the BL Lac 
objects, our results are very similar to Ledden and O'Dell's (1985) results 
because of the similar database. We note, however, substantial differences 
among different studies of optically selected quasars, ranging from <a^^>-l .37 
± 0.10 (Marshall et al. 1983) to <q ox >— 1.65 ± 0.03 (this study). Since 
Marshall et al. treat lower luminosity objects, this may cause the difference 
as we can see in an other study (Zamorani et al . (1981) find for radio quiet 
quasars that ct^-1 . 37+0 . 05/-0 . 08 for log(i o )<31.4 but a =1 . 62+0 . 08/-0 . 11 for 
log(i o )>31 .4) . Other possible causes for discrepances are the use of 
different optical magnitudes (i.e. not those of Hewitt and Burbidge 1980), 
frequences at each band (Owen, Helfand, and Spangler 1981, Cruz -Gonzales and 
Huchra 1984), selection criteria (Zamorani et al . 1981), correction factors 
for absorption and Mgll line (Margon, Downes, and Chanan 1985), and different 
spectral indices assumed for extrapolations. 
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c) Luminosity Ratio Dependence on Optical Luminosity and Redshift 


The dependence of a on optical luminosity and/or redshift for the 
optically selected quasars is often discussed (Reichert et al. 1982, Avni and 
Tananbaum 1982, Zamorani 1982, Maccacaro and Gioia 1983, Zamorani 1984, Kriss 
and Canizares 1985, Avni and Tananbaum 1986). Avni and Tananbaum (1982) were 
the first to obtain a relation among them using survival analysis (the 
" detection- and-bounds" linear regression method, see Avni and Tananbaum 1986). 
Their relation is expressed as 

a ox -(-° . 0±0 . 3) (r-0 . 5 )+(0 . 12±0 . 06) [ log(i o )-30 . 5 ] +1 . 50 , 
where r - z/(l+z) . They mention that explicit dependence of a ^ on the 
optical luminosity is predominant, but the joint dependence of on both 
variables is possible. This result is confirmed in a more recent paper (Avni 
and Tananbaum 1986). 

Zamorani (1982) shows a similar relation for a combined sample of 

optically selected quasars and Seyfert galaxies . Ignoring the dependence on 

the redshift, Zamorani (1984) finds 

a -0.129 log (2 )-2 .427 . 
ox & o 

In another subsample, Tanambaum, Wandle , and Zamorani (1983) find a similar 
relation for radio selected quasars (the 3CR sample) , 

log(i x )-27.63+(0.47±0.15) [log(2 o )-31 . 20]+(0 . 14±0 . 12) [log(i r )-34. 78] 

-(0 .45±0 .78) [log(l+z)-0 . 26] . 

Because of the weak dependence on the radio luminosity and the redshift, they 
rewrite this relation as 

a =0.20 iog(i )— 4.98. 
ox ° o 

We look for analogous relations in our data sets; the results are shown 

in Tables 9a and 9b. In Table 9a, we show the relations between a and r and 

ox 
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between a and log (i ) . Assumed regression forms are a =*a+br and 
ox ° o ox 

a =a+blog(i ). We find the log(i ) - a slope to be 0.11 ±0.02 for radio 
selected and optically selected quasars. Only the optically selected quasars 
show high significance levels for both the redshifts and the optical 
luminosity regressions. These relations, however, might be artificial, since 
for the optically selected quasars, the optical luminosity is biased due to 
optical magnitude limited survey . The direction of evolution of the optically 
selected quasars agrees with the results by Ku, Helfand, and Lucy (1980) and 
Zamorani et al . (1981) (see Table 1). 

We also compute the three dimensional regressions for a , T , and 

log(i o ) , using the regression form of Avni and Tananbaum (1982). The results 

are shown in Table 9b. The second column shows Cox probabilities. The first 

value is a joint probability that no correlation exists between a and both 

ox 

redshift and log , the second value is the probability for the redshift 
alone, and the third is for optical luminosity density alone. The second and 
third values are determined assuming the ratio of the slope coefficient and 
the error is distributed as a Gaussian. The third column shows the regression 
results. All subsamples, except the X-ray selected BL Lac objects (P-23%) , 
show highly significant joint probabilities (P < 0.01%). Only the optically 
selected quasar sample does not show significant correlations for the 
individual variables (P(z-)=53% and P(i Q )=42%) , even though this subsample is 
the only one which shows high significance levels for correlation between a ^ 
and r, and a and log(i Q ) . Although this result does not confirm Avni and 
Tananbaum' s result which shows that a QX is positively correlated with log(2°) , 
we find similar relations (high joint and individual probabilities) in other 
subsamples . 
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The interpretation for the radio selected quasars, X-ray selected AGNs , 
and BL Lac objects is that objects at higher redshifts have higher 
ratios, and those with higher optical luminosities have lower H x /% 0 ratios. 

For example, a typical radio selected quasar at z-2 will have a 2. x /& 0 ratio 
twice that of a similar quasar at z-0. At a given redshift, a typical radio 
selected quasar with log(i o )-32 will have a 2 ratio half that of a similar 
quasar with log(i o )-33. The slope coefficients in Table 9b can be used to 
give analogous results for other subsamples . 

We also examine the dependence of a on an X-ray luminosity and a 
redshift. The results are in Table 9c. Because of the presence of censoring 
in as well as a , we can obtain results only for the X-ray selected AGNs 
and the BL Lac objects. Although the X-ray selected AGNs show a high 
significance level for both the variables, the BL Lac objects do not. 

Comparing these three results, we find some inconsistencies. The radio 
quasars, X-ray AGNs, and BL Lac objects show weak positive correlations 
between and r in Table 9a, but strong negative correlations in Table 9b. 
Table 9c shows another problem. In the relations among a , r, and log(i ), 
the direction of the dependence on r is opposite to that of r in the relations 
among <* ox > r, and logCi^) . If there is a real dependence of a on r, the 
direction should be the same. These inconsistencies suggest that either a 

ox 

does not truly depend on r , or that the evolution is different for the various 
wavebands and sub samples . 

d) Two Component Model for X-ray Emission of the Radio Selected Quasars 

Although an important emission mechanism of the radio selected quasars is 
thought to be the synchrotron radiation, the slope of the log(i^) - log(i ) is 
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much shallower (Z —Jp'^ ^'^) t h at Q f gg, L ac objects (Z _^0.77+0.08) ^ 

whose emission mechanism is almost certainly the synchrotron. A close 
examination of the plot (Fig. la) shows that the distribution of data points 
does not follow a straight line. This has been be interpreted as evidence for 
two different types of the radio selected quasars (Owen, Helfand, and 
Spangler 1981, Owen and Puschell 1982, and Zamorani 1984). 

Zamorani (1984) has suggested that there are two different X-ray emission 

mechanisms; for example, a synchrotron component and a thermal component. For 

the relation between Z^ and Z (see Fig. la), the emission mechanism of the 

steeper component at high luminosities would be mainly non-thermal, and that 

of the flatter component at lower luminosities would be thermal or a 

combination of the non-thermal and thermal emission. Since if the X-ray 

emission is a thermal origin, the X-ray luminosity is expected to be 

independent from the radio luminosity. This explanation is supported by the 

partial correlation coefficients studied in the § Va, since the partial rank 

coefficients between the radio luminosity and the X-ray luminosity and that 

between the optical luminosity and the X-ray luminosity are equally strong. 

The X-ray emission related to the radio luminosity may be a synchrotron or SSC 

radiation because of the similarities to the BL Lac objects, and the emission 

related to the optical luminosity may be a thermal radiation because of the 

similarities to the optically selected quasars. Zamorani adopts a three 

dimensional regression model with a form of 

b b 

Z oc 2 ° + Z r . 
x o r 

He finds that for the flat spectrum radio selected quasars, 5^=0. 75 and 
0.95, and for the steep spectral radio selected quasars, 5^=0. 63 and 
b^-0.75. Since the survival analysis cannot treat non-linear problems, we 
adopt Zamorani' s b-values. Then, using our data set and the application of 
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the BHK method, we find his equations can be rewritten as 

log(i x )-log(/- 75 +5xl0" 1 V- 95 )+21.8, 
log ( -log ( i° ' ■' 6 3 +3xl0~ 7 Jt° r ■ 1 5 ) +2 6 . 0 , 
respectively. The model seems to fit well; however, we found that the 
computation of b -values heavily depends on a few points and hence the value is 
unstable . 

To judge whether this two component model is needed to fit the data, we 

compare these models to the other simple models. One is a straight line 

(log(i x )-29 . 0+0 . 481og(i^_) ) and another consists of two straight lines 

(log(i )-19 . 0+0 . 771og(i ) for log(i )>33.77 and log(je )-45.0 for 
x i* r x 

l°g(i r )<33 . 77) . The last model is made by assuming that the steeper 

components have the same slope coefficient as the BL Lac objects and the flat 

components have a zero slope coefficient. The other coefficients are then 

2 

found using the BHK method. Since we do not have either a x -test or an other 
goodness -of- fit test to compare models for censored data, we need to use a 
non-standard method. The Kaplan-Meier estimation of distribution of residuals 
found by subtracting these models from the data is examined. The 25th and 
75th percentiles of the residuals express the dispersion of the data about the 
model. One problem is that these dispersions cannot be translated to 
probabilities, and hence the results are only qualitative. 

We find dispersions of 0.65, 0.63, 0.51, and 0.44 about the straight line 
model, the two straight line model, the flat radio spectra model, and the 
steep radio spectra model respectively. Zamorani's (1984) models are thus 
better than these other models of the X-ray emission for radio selected 
quasars . 


Additional support for these composite 
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unresolved radio quasars . The subsample of the resolved quasars (with 
arcsecond resolution) has a shallower slope (i compared to the 

subsample of the unresolved quasars (i 71-0.07 £ or unresolved quasars 

with £*^<0.3; see Fig. 3 in Paper II). These results are confirmed in our 
enlarged data set using survival analysis. 

If the X-ray emission mechanism of the resolved radio selected quasars is 
dominantly thermal , and that of the unresolved radio selected quasars is non- 
thermal, then we can expect a flat slope for the resolved, and a steep slope 
for the unresolved radio selected quasars . 

e) BL Lac Objects 

In our data sets, the radio selected BL Lac objects have a distinct 
position. Since all the objects are detected in all three frequencies, 
statistical results are free from the selection effect due to the flux limited 
observations. Also because of their nonthermal nature (supported by short 
variabilities, high polarizations, and our partial correlation results), they 
can be used as a standard to which other subsamples are compared. 

Since the BL Lac objects have no upper limits, the partial linear 
correlation probabilities in Table 7 can be computed to find which relations 
are significant. The partial linear correlation probabilities are 
P(rx, o)<0 . 01% , P(ox,r)=2%, and P(ro,x)«20% , where, for instance, P(ro,x) is 
the correlation probability between the radio and the optical luminosity for a 
fixed X-ray luminosity. The results show that the radio/X-ray relation is 
most significant, and the optical/X-ray relation is moderately significant, 
but the optical/radio relation is not significant. 
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Some researchers have pointed out the similarities between the BL Lac 
objects and the radio selected quasars. Using two sample tests, we compare 
these subsamples. First, we use all data in both the samples. The results 
are shown in Table 10a. The optical and the X-ray luminosities of the radio 
selected quasars are higher than those of the BL Lac objects, which is 
probably a consequence of the fact that the radio selected quasars have a much 
wider distribution in the redshift than the BL Lac objects. If the redshift 
range for both the subsamples is restricted to 0.8 and 1.7, better agreements 
in the mean luminosities are found (Table 10b). There are, however, some 
problems. There are only eight BL Lac objects in this restricted sample, and 
they do not show any internal radio, optical, or X-ray correlations. The 
radio selected quasars, in contrast, give the slope coefficients 0.34+0.10/- 
0.08 for log(f )-log(i ) and 0.45 ± 0.11 for log(i )-log(i ), though no 
significant relation exists between log(i o ) and log(i^) . These results can be 
interpreted in two ways. One interpretation is that there are no significant 
differences in the continuum spectra of radio selected quasars and BL Lac 
objects, and another is that the data are too fragmentary to give firm 
conclusions. More BL Lac objects are needed. 

The unresolved radio selected quasars are also compared with the BL Lac 
objects. The samples with restricted redshift range show strong similarities. 
Therefore the BL Lac objects may be special cases of the radio selected 
quasars . 

The distinctive difference between the radio selected and the X-ray 

selected BL Lac objects are often noted (Ledden and O'Dell 1985, Stocke et al . 

1985). For example, the mean spectral indices are: <c*^ o >= 0.62 ± 0.02, 

<a >=1.46 ± 0.03 for the radio selected BL Lac objects, and <a >=0.37 ± 
ox J ro 

0.02, <a >=1.11 ± 0.04 for the X-ray selected BL Lac obiects. Ledden and 

OX " ~ 
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O'Dell (1985) suggest that the main emission mechanisms of the X-ray BL Lac 
objects is the synchrotron radiation, and the radio selected BL Lac objects 
have extra mechanisms, such as beaming. If we compare the averaged 
luminosities of the radio selected and the X-ray selected BL Lac objects, we 
find that the radio selected BL Lac objects are 100 times more luminous than 
the X-ray selected BL Lac objects in the radio band, 7 times more luminous in 
the optical band, and nearly same in the X-ray band. This may support the 
beaming model . 

But an alternative possible explanation is a selection effect. In the 
diagram of a ro / a ox relation, we see that the radio selected and the X-ray 
selected BL Lac objects mark the lower and upper bounds of the radio selected 
quasars . The difference between these two groups may be due to two extreme 
cases of sources with the same emission mechanism as in the radio selected 
quasars. A similar offset in the <a rQ > is seen between the radio selected 
quasars and the optically selected quasars and the X-ray selected AGNs which 
can be attributed largely to the selection methods used in their discovery. 

If so, we may find "optically selected" BL Lac objects somewhere between these 
two groups. Only one BL Lac object was possibly found optically (ZWI 186), 
and it is located among the X-ray selected BL Lac objects. Optical surveys 
for high polarized objects (e.g. Borra and Corriveau 1984) have been generally 
unsuccessful. 

f) Comparison of the Optically Selected Quasars and the X-ray Selected AGNs 

It is often mentioned that the optically selected quasars and the X-ray 
selected AGNs have similar natures, and the latter are treated as a lower 
luminosity sequence of the former (e.g. Maccacaro et al. 1984, Kriss and 
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Canizares 1985). In our data sets, both subsamples show a strong correlation 
between the optical and the X-ray luminosities but do not show the other 
relations (see § Va) . This initially suggests that the X-ray and optical 
emission mechanisms of these subsamples are similar. In the plot of log (i ) 
and log(X x ) , however, the slope coefficient of the X-ray selected AGNs are 
significantly steeper (0.87 ± 0.05) than that of the optically selected 
quasars (0.70 ± 0.06). Also the indices of the X-ray selected AGNs tend 

to be shallower ( <a ox > ”1-65 ± 0.03) than those of the optically selected 
quasars (<a^>^ 1.35 ± 0.02). The average radio and optical luminosities of 
the optical selected quasars are about 20 times brighter than those of the X- 
ray selected AGNs, but the average X-ray luminosity is nearly the same. If 
the emission mechanism of these two subsamples were the same, the difference 
of the brightness of each band should be approximately same. One explanation 
of this difference is suggested by Kriss and Canizares (1985). They show that 
the high redshift objects are much "redder", since the optical band shifts to 
shorter wavelength which is strongly affected by reddening, while the X-ray 
band is not affected much. Hence, if we compare the optically selected 
quasars which have, on the average, a higher redshift to the X-ray selected 
AGNs which have, on the average, a lower redshift, then the optically selected 
quasars show more absorption. 


(g) Groping towards the Physics of AGN Continua 


Having investigated the relations between the continua of various types 
of AGNs, we should like to know the relevance of the various models of 
physical processes of continuum emission to the results of our statistical 
analyses. There are fundamentally two different mechanisms : a thermal 
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radiation from an accretion disk, and a non- thermal radiation from the 

vicinity of the central engine or jets. If the thermal emission from the 

accretion disk is the main mechanism (Tucker 1983, Schlosman, Shaham, and 

Shaviv 1984) , X-ray and optical luminosities may show a correlation but the 

radio luminosity is likely to be independent. According to Tucker (1983), the 

emission from an optically thick accretion disk can generate the relation 

i oc % a with a—0.5 to 0.8. Most of our results are consistent with this 
x o 

prediction (a-0.7, except for BL Lac objects where a - 0.9), though the 
thermal model does not explain the radio/X-ray correlation seen in most 
samples. If the entire continuum is due to synchrotron emission, the X-ray, 
optical, and radio luminosities should be well correlated each other with 
a o ^a ro . All subsamples agree with these conditions. If the synchrotron 
self-Compton (SSC) mechanism is important, a strong correlation between the 
radio and X-ray is expected with a possible correlation with the optical 
luminosity through the synchrotron emission. If beaming due to a relativistic 
jet exists, it would lead to correlation between all beamed (presumebly 
nonthermal) components. Other mechanisms, such as Compton scattering of 
blackbody or cyclotron radiation could also be responsible for the power low 
spectrum in the optical to X-ray bands. 

These various models clearly do not make predictions which can be 
uniquely distingushed by the Z J 2 database studied here. Most models are 
not sufficiently developed to predict how radio, optical and X-ray 
luminosities should scale. Nontheless , we can attempt to reach some crude 
conclusions. The fact that 2^ is correlated with both 2 q and 2^ in virtually 
all subsamples of AGNs (Table 5a) is evidence against a thermal model for the 
continuum spectrum unless, for example, there is some indepedent scaling 
between the size of the thermal accretion disk and the strength of the non- 
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thermal jets. The correlation between all three bands is fully consistent 
with a simple or beamed synchrotron model, though models must account for fact 
that the scaling between bands is not quite linear (Table 5b) . There are no 
indications that £ and Z^ are correlated with £ decoupled, as might occur in 
an SSC model, but this cannot be conclusive evidence against SSC as the 
optical band could be dominated by non- thermal continuum (e.g. Konigl 1981). 

Although it is risky to pursue more elaborate models when adjudication 
between the simplest ones is difficult, we find the two component AGN 
continuum model discussed in § Vd is attractive. Here all AGNs have the 
thermal emission from the accretion disk and the non- thermal emission from the 
jets. The differences between subsamples may be due to the differences 
between accretion modes (Blandford 1984) . If a radiation torus around a black 
hole is radiating at just over the Eddington limit, the emission is dominated 
by the thermal radiation, since there are insufficient relativistic electrons 
to power a synchrotron continuum. The radio loud quasars accrete at higher 
rate so that they produce the jets populated by relativistic electrons, but 
the synchrotron emission (and possibly SSC) need not dominate and hence we see 
both the thermal disk and non- thermal jet radiations. The BL Lac objects may 
be an extreme case with intrinsically luminouse jets and a faint optically 
thin accretion disk, or they may possess ordinary jets that happen to be 
pointed to us so that the synchrotron emission is extremely enhanced and 
dominates the total luminosity. 

There are a few complications in any of these interpretations. One 
concern is possible evolution effects. If the continuum emission in various 
bands are intrinsically related, then correlations should appear even if the 
subsamples are divided into narrow redshift bins within which no evolution 
could occur. We find most of the Z j£ Jz ^ correlations are present in 
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specified bins , but appear weaker than the correlations seen in the entire 
subsamples. For example, the correlation probabilities between the radio and 
the X-ray luminosities of the radio selected quasars are 5% for 0.0 < z < 

0.5, 0.01% for 0.5 < z < 1.0, 0.03% for 1.0 < z < 1.5, 93% for 1.5 < z < 2.0, 
17% for 2.0 < z < 2.5, and 0.6% for 2.5 < z. The low correlation 
probabilities in these subsamples are partly due to the reduced size of the 
data sets in each bin, and partly by the narrow range of luminosities in each 
redshift bin. The existence of the correlations within narrow redshift ranges 
gives some confidence, however, that the correlations are not entierly due to 
cosmological luminosity evolution. 

Another problem is inappropriately defined samples. We subdivided our 
samples by selection criteria such "X-ray selected" or "optically selected" 
objects. This method may introduce some mixing of intrinsically different 
samples . For example , the optically selected quasars clearly contain a few 
"radio loud" quasars. Although this mixing might lead to some misleading 
statistical results, comparison between the "radio selected" (which are chosen 
from their initial discovery in radio surveys) and the "radio loud" subsamples 
(which are chosen from the entire sample with the condition a >0.3) shows 
that there is no significant difference. 

We thus find that, although a large number of data were collected, it 

proves difficult to specify a physical model. This is partly because most 

theoretical studies do not show tracks in a Ja or Z / £ / H plots. 

ox? ro r ' o x 

Since these kinds of plots are now widely produced observationally , we 
encourage theorists to make such predictions. It is also desirable to have 
deeper surveys in all bands so that we can examine samples with wide 
luminosity ranges within specific redshift ranges. These surveys may clarify 
the evolution effect on the emission mechanisms. 
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VI. SUMMARY 


Using a large database, we investigated statistical properties of AGNs 
continuum levels in the radio, optical, and X-ray bands. For the radio 
luminosity of AGNs, we used the core luminosity to discount effects from radio 
lobes and jets. The statistical methods called survival analysis were used to 
show the statistical relations despite the upper limits in flux limited data 
sets. Our main results are as follows: 

1. For the optically selected quasars and the X-ray selected AGNs, Z^ is 

correlated with 1 , but not with Z . Also Z does not correlate with Z . This 
x r r x 

suggests that the main emission mechanism of these subsamples is thermal 
emission. For the BL Lac objects, ^ r /^ x relation is most significant, and 
this suggests that the main emission mechanism of the BL Lac objects is non- 
thermal. The radio selected quasars show high significance levels for both 
Z r /^ x and ^ Q /^ x relations. Hence they may have both the mechanisms, which is 
further supported by model fitting suggesting that the radio selected quasars 
have two different X-ray emission mechanisms. 

2. The BL Lac objects have similar emission mechanisms as the radio selected 
quasars, at least in the limited redshift range overlapping both samples. The 
radio selected BL Lac objects are perhaps special cases of the unresolved 
radio selected quasars. The difference between the radio selected and the X- 
ray selected BL Lac objects may be due either to beaming effects or selection 
effects . 

3. Some compact radio sources with faint optical counterparts show that 
<a^> > <c* ox >. This suggests the possibility that SSC emission may be present 
in the optical to X-ray bands. 

4. The spectral index between the optical and X-ray luminosities depends on 
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the optical luminosity. The optical luminosity increases faster than the X- 
ray luminosity. 
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Table Captions 


Notes to Table 2a: 

1. First detected by Condon and Dressel (1978). We find a faint radio 
lobe (3 raJy) 6" at P.A. 104 ° from the quasar. 

2. A 6xRMS detection within 1" of the optical position. 

3. A 5 . 5xRMS detection within 1" of the optical position. Mrk 205 was 
also detected by Sulentic (1986) at a level of 1.48 ± 0.25 mJy 

at 5 GHz . 

4. First detected by Sramek and Weedman (1980). Our improved position is 
12h58m59.4s, 34°16'38" . 

5. A 7xRMS detection within 1" of the optical position. 

6. First detected by Sramek and Weedman (1980). Our improved position is 
16h04m53.4s, 29°03'21". 

Definitions for Table 3: 


Optical 

If a radio AGN has an optical "empty field", V>20.0 is used as the optical 
upper limit. 

X-ray 

If only the X-ray flux (or flux density or count number) is available, the 
luminosity is computed according to the descriptions below for the given 
reference : 

* LI X-ray list gives X-ray data in the flux density at 1 keV. The conversion 
to the flux (0.5 -4. 5 keV) is 
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S x (0.5-4.5) [10 13 erg/sec/cm 2 ]= 0.068 f x (lkeV) nJy 

* 01 X-ray list gives the X-ray data in the flux (0.15-3.5 keV) . The 
conversion is 

S (0 . 5-4 . 5)=0 . 95 S (0.15-3.5) 

X x 

* B3 X-ray list gives the X-ray (0. 5-3.0 keV flux). The conversion is 

S (0 . 5-4 . 5)=1 . 38 S (0. 5-3.0). 

X X 

* G2 X-ray list gives the X-ray (0.3-3. 5 keV flux). The conversion is 

S (0.5-4. 5)=1. 07 S (0.3-3. 5). 

•X x 

* W2 X-ray list gives the X-ray count rate. The conversion is 

S^-2. 88x10 13 X 2 erg/ sec/ cm' . 

* The X-ray flux is computed from Einstein 1PC photon counts (cts/sec) by 

1.0IPC(cts/sec)-3.0xl0 ^ erg/sec/cm' 

20 —2 _n s 

assuming N(H)=* 3x10 cm and S —( freq ) ' . 

* The X-ray flux density at 2KeV is computed by 

2 

S (2KeV)-l.t*7xlO~ 18 f er f/f ec /^ m E 
* * £ 1°' 5 _ £2 °- 5 

o 

where El and E2 are the band limits in KeV and f is in erg/ sec/ cm . 

Notes to Table 8: 

1. Data have too many upper limits and the result is obtained from 
a limited range. 
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Table 1 : Spectral Indices from Recent Literatures 

Study Spectral Indices Descriptions 

Tananbaum et al . a -1. 2710.07 Mixed QSOs 

ox x 

(1979) 

Ku et al . (1980) a —1 .46±0. 02 Total sample 

a — 1. 3810.03 Radio selected QSOs 

ox x 

a — 1 . 5210 . 03 Optically selected QSOs 

a “1.4110.03 X-ray selected AGNs 

ox 

a ox ’~ 1 ■ 3610 . 04 Radio QSOs with low redshift (z<1.0) 

® ox ”l. 4-010.03 Radio QSOs with high redshift(z>l . 0) 

« ”1 • 3610. 04 Opt. QSOs with low redshift (z<1.0) 

a ox ”1.6510.04 Opt. QSOs with high redshift (z>1.0) 

a “1.2510.05 OW 

ox 

a —1.3110.05 BL Lacs 

ox 

Zamorani et al. a -1.2710.03 Radio loud 

ox 

(1981) a —1.46+0.05/-0.07 Radio quiet 

ox 

« ox -1.35+0.05/-0.08 Radio quiet with low redshift (z<1.0) 
c* ox ”l . 62+0 . 08/-0 . 16 Radio quiet with high redshift (z>1.0) 
a ox “l . 37+0 . 05/-0 . 08 Radio quiet with log(j? )<31.4 
a ox =l . 62+0 . 08/-0 . 11 Radio quiet with log(,0 )>31.4 


Owen et al. 

(1981) 

a =1.0210.05 
nix 

mm selected AGNs 



a “1.2110.19 
ox 

mm selected AGNs 

Stocke et al 

. (1983) 

a =1.310.2 
ox 

X-ray selected AGNs 

Marshall et 

al. 

a =1.3710.10 
ox 

Optically selected quasars 


(1983) 
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Cruz -Gonzales and 


Huchra (1984) 
Margon and Chanan 
(1985) 


a 0 . 59±0 . 10 
ri 

BL Lac 

: radio 

- infrared 

a -0 . 94±0 . 09 
rx 

BL Lac 

: radio 

- X-ray 

a -1 . 26±0 . 03 
ox 

X-ray 

selected 

quasars 


Ledden and O'Dell 

a -0 . 63±0 . 12 
ro 

BL Lacs 

(1985) 

a -0 . 89±0 . 06 
rx 

BL Lacs 


a -1. 4010.17 
ox 

BL Lacs 


a =0 . 74±0 . 10 
ro 

HPQs 


a -0.9210.06 
rx 

HPQs 


a -1.3010.13 
ox 

HPQs 


a -0.6710.12 
ro 

Blazars 


a -1.3610.17 
ox 

Blazars 
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Table 2a : VIA Observations of Optically Selected Quasars 


Obj ect 

Name 

Sg(mJy) 

Note 

Obj ect 

Name 

S 6 (mJy) 

Note 

0137-010 

NAB 

<0.8 


1045+128e 

nrNGC3384 

<1.1 


0143-015 

MC5 366 

<0.8 


1045+128f 

nrNGC3384 

<0.9 


0143-010 

MC5 368 

<0.8 


1045+128g 

nrNGC3384 

<0.9 


0146+017 

MC5 141 

<1.2 


1045+128h 

nrNGC3384 

<1.0 


0207-378 


<1.3 


1202+281 

GQ Com 

1.1 

2 

0241+011a 

nrNGC1073 

<2.0 


1219+755 

Mrk 205 

0.9 

3 

0241+011b 

nrNGC1073 

<2.0 


1246-057 


<0.9 


0241+Ollc 

nrNGC1073 

40.5 

1 

1258+286 

W 61972 

<0.9 


0242-410 


<1.4 


1258+342 

KP 33 

25.1 

4 

0849+154 

LB 8796 

<0.9 


1300-243 


1.3 

5 

0854+194 

LB 8948 

1.8 


1334+286 

RS 23 

<0.8 


0855+188 

LB 8991 

<0.9 


1346-036 


<1.0 


0856+186 

LB 9010 

<1.0 


1604+290 

KP 63 

4.0 

6 

0856+189 

LB 9029 

<0.8 


1606+288 

KP 64 

<0.9 


1045+128a 

nrNGC3384 

<1.0 


1606+289 

KP 67 

<0.9 


1045+128b 

nrNGC3384 

<0.9 


1720+246 

V396 Her 

31.0 


1045+128c 

nrNGC3384 

<0.9 


1803+676 


<0.8 


1045+128d 

nrNGC3384 

<1.0 


2225-055 

PHL 5200 

<0.8 
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Table 2b : VIA Observations of X-ray Selected AGNs . 


Obj ect 

Name 

Sg(mJy) 

Object 

Name 

Sg (mJy) 

0057+311 

IE 

<0.7 

1205+465 

IE 

<0.7 

0112+325 

IE 

<0.8 

1228+164 

IE 

<0.8 

0225+312 

IE 

<0.6 

1304+341 

IE 

<0.7 

0244+192 

IE 

<0.6 

1352+1820 

IE 

<0.7 

0357+104 

IE 

<0.7 

1352+1828 

IE 

<1.0 

0745+554 

IE 

<0.7 

1357-022 

IE 

<0.7 

0754+392 

IE 

3.6 

1401+095 

IE 

18.8 

0906+425 

IE 

<0.8 

1529+050 

IE 

3.8 

1008+345 

IE 

<0.8 

1530+151 

IE 

<0.9 

1011+032 

IE 

<0.7 

1602+241 

IE 

<1.2 

1031+582 

IE 

<0.6 

1747+683 

IE 

<0.6 

1139+104 

IE 

<1.1 

2251-175 

IE 

<0.7 
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(c) X-ray Selected Active Galactic Nuclei 
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Table 6a Correlations 


Sample 


Total 

Censored 

Significnat level 

Variables no. 

obj ects 

of correlations Plot 

objects 


no . 



X 

y 


X 

y 

both 

Cox (%) 

BHK (%) 


Radio Selected 

log(J? r 

)-log(f x ) 

156 

50 

7 

4 


<0.01 

la 

Quasars 

lo g(i o 

)-log(i x ) 

156 

0 

11 

0 

<0.01 

<0.01 

2a 


lQ g(i 0 

)-log(i r ) 

156 

0 

54 

0 

<0.01 

<0.01 

3a 


a 

ro 

a 

ox 

156 

50 

7 

4 

<0.01 

<0.01 

4a 

Optically 

log (i r 

)-log(i x ) 

103 

18 

6 

68 

. . . 

<0.01 

lb 

Selected 

10 g(i 0 

)-log(i x ) 

103 

0 

74 

0 

<0.01 

<0.01 

2b 

Quasars 

lo «< J o 

)-log(i x ) 

103 

0 

86 

0 

0.02 

0.2 

3b 


a 

ro 

a 

ox 

103 

18 

6 

68 


<0.01 

4b 

X-ray Selected 

i°g(^ r 

,)-log (i x ) 

122 

103 

0 

0 


<0.01 

lc 

AGNs 

io g ( i o 

)-log(i x ) 

122 

0 

0 

0 

<0.01 

<0.01 

2c 


log(i Q 

)-log(i r ) 

122 

0 

103 

0 

<0.01 

<0.01 

3c 


a 

ro 

a 

ox 

122 

103 

0 

0 


52 

4c 

Radio Selected 

log(^ r 

.)-io g (i x ) 

24 

0 

0 

0 

<0.01 

<0.01 

Id 

BL Lac Objects 

10 &(* 0 

,)-log(i x ) 

24 

0 

0 

0 

<0.01 

<0.01 

2d 


log(^ 0 

,)-log(i r ) 

24 

0 

0 

0 

<0.01 

<0.01 

3d 


a 

ro 

a 

ox 

49 

0 

2 

0 


<0.01 

4d 


50 



X-ray Selected 

log(i r )-log(i x ) 

10 

0 

0 

0 

2 

3 

Id 

BL Lac Objects 

log(i )-log(i ) 

O -A 

10 

0 

0 

0 

0.03 

0.3 

2d 


log(X o )-log(i r ) 

10 

0 

0 

0 

0.7 

0.7 

3d 


a - a 

ro ox 

16 

0 

0 

0 

4 

2 

4d 

Compact Radio 

a - a 

ro ox 

19 

0 

1 

8 


0.03 

4e 

Sources 









Flat Spectral 

log(J r )-log(J x ) 

66 

0 

2 

0 

<0.01 

<0.01 


Radio Quasars 

log(i o )-log(i x ) 

66 

0 

2 

0 

<0.01 

<0.01 



log(i o )-log(ip 

66 

0 

0 

0 

<0.01 

<0.01 



a - a 

ro ox 

66 

0 

2 

0 

<0.01 

<0.01 

. . 
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Table 6b Linear Regressions 


Sample 


Variables 


Total Censored 

no. objects 

objects 


x y 


x y both 


Intercept Coeff. 
Slope Coeff. 
Standard Deviation 
EM B-J 


Radio Selected log(i^) -logCi^) 156 50 

Quasars 


4 29.0±... 

0.4810.06 


Optically 

Selected 

Quasars 


log(i o )-log(i x ) 

156 

0 

11 

0 

23.911.6 

+i 

00 

CO 

CM 






0.7010.05 

0.7010.05 






0.47 

0.45 

log(i o )-log(i r ) 

156 

0 

54 

0 

5.6213.73 

6.051. . . 






0.9110.12 

0.8910.12 






1.05 

0.96 

a - a 

ro ox 

156 

50 

7 

4 

-0.7510 

.12 

i °g(- ? r )- l °g(^ x ) 

103 

18 

6 

68 

28.21 







0.51+0 . 15/-0 . 11 

log(i )-log(i ) 

O A 

103 

0 

74 

0 

23.1+1.9 

23.11 ... 


0.7010.06 0.7010.07 

0.50 0.48 


52 


t 


log(i )-log(i ) 103 0 86 0 0.3210.65 -0.681 ... 

0.9710.21 1.0110.27 

1.39 1.25 

a - a 103 18 4 68 

ro ox 

-0 . 60+0 . 04/-0 . 05 


X-ray Selected log(i^) -log(i^) 122 103 0 0 

AGNs 


0.46+0 . 08/-0 . 06 


log(i )-log(i ) 122 000 18.211.4 

0.8710.05 

0.43 

log (£ )-log(i ) 122 0 103 0 -34.5112.1 

o r 

2.0810.40 

2.42 

a -a 122 103 0 0 

ro ox 


Radio Selected log(i^) -logCi^) 24 000 

BL Lac Objects 

log(je o )-log(i x ) 24 0 0 0 


18.712.6 

0.7710.08 

0.50 

17.813.3 

0.8910.11 

0.59 


53 


0 


0 


0 


3 . 2314 . 13 


log(i Q )-log(i r ) 24 

1 . 00±0 . 14 
0.74 

“ro ' “ox 49 020 1 . 95±0 . 11 1.95±... 

-0 . 93±0 . 18 -0 . 93±0 . 19 

0.18 0.18 


X-ray Selected log(i^) -log(^) 10 0 0 0 

BL Lac Objects 


log(i o )-log(i x ) 10 0 0 0 


log(i o )-log(i r ) 10 0 0 0 


a - a 
ro ox 


16 0 0 0 


14 . 4±0 . 92 
0 . 96±0 . 29 
0.34 
11 . 4 ± 6 . 9 
1 . 13±0 . 23 
0.27 

14 . 7±7 . 9 
0 . 56±0 . 27 
0.31 

1 . 59±0 . 15 
-1.40±0.41 
0.11 


Compact Radio a -a 

ro ox 


19 0 1 8 


Sourcs 


-0 . 97+0 . 23/-0 . 16 


Flat Spectral log(i^) -log(i x > 66 
Radio Quasars 


log(J! o )-log(i x ) 66 


0 

2 

0 

22 . 6±2 . 7 

22. 6±. . . 




0 . 67±0 . 06 

0.6710.08 




0.43 

0.42 

0 

• 2 

0 

24 . 2±2 . 0 

24.01. . . 




0 . 69±0 . 06 

0.7010.06 




0.42 

0.39 


54 



66 0 


0 


0 


11 . 5 ± 2 . 7 


log(i o )-log(i r ) 

0 . 7410 . 01 
0.58 

a ro ' “ox 66 0 20 1.9110.11 1.981... 

-0.9210.17 -0.9110.15 
0.16 0.15 
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Table 7 : Partial Rank Correlation Coefficient 


Sample Correlation 

Fixed variable 


Radio QSOs 
Optical QSOs 
X-ray QSOs 
Radio BL Lacs 
X-ray BL Lacs 


) -log(i ) log(lo)-log(i ) log(i )-log(i ) 

-f 2Z X. O L 

log(i fl ) log(^p log(i x ) 


0.31 

0.45 

0.13 

0.05 

0.16 

0.08 

0.07 

0.64 

0.07 

0.56 

0.45 

0.34 

0.09 

0.60 

0.48 


Table 8: Kaplan-Meier Mean and Standard Deviation 


of a and a 
ro ox 


Data 

<o > 
ro 

<Q > 

OX 

Radio QSOs 

0 . 6410 . 01 

1.43+0.02/-0.03 

Optical QSOs 

C0.49t0.02) 1 

1.6510.03 

X-ray AGNs 

0.48+0. 05/-0. 06 

1.3910.02 

Radio BL Lacs 

0 . 62±0 . 02 

1.4610.03 

X-ray BL Lacs 

0 . 37±0 . 02 

1.1110.04 

Compact Radio Sources 

0 . 93±0 . 03 

1.1210.03 

Flat Spectrum Radio 

0.6310.02 

1 . 39+0 . 03/-0 . 04 
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Table 9a 


: Correlations between a and r or lo z(£ ) 

ox b o 

t - a log(i ) - a 

ox ° o ox 

Sample Cox prob. Intercept Cox prob. Intercept 

censored no. Slope Slope 

(%) Stand. Dev. (%) Stand. Dev. 

total no. 


Radio QSOs 
(11/156) 

9.1 

1.27±0.05 
0 . 21±0. 10 
0.20 

<0.01 

-2.2110.62 

0.1110.02 

0.19 

Optical QSOs 
(73/102) 

<0.01 

1 . 35±0. 07 
0 . 54±0. 13 
0.21 

<0.01 

-1.9410.72 

0.1110.02 

0.19 

X-ray AGNs 
(0/123) 

37 

1.3810.03 

-0.1310.09 

0.17 

0.13 

-0.1010.54 

0.0510.02 

0.16 

Radio 
BL Lacs 
(0/24) 

34 

1.4710.08 

-0.2510.22 

0.22 

29 

0.0511.25 
0 . 0410 . 04 
0.22 

X-ray 
BL Lacs 

(0/10) 

9.7 

1.0710.05 

-0.3910.32 

0.10 

85 

0.0012.72 

0.0310.09 

0.11 
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Table 9b : Correlations among a , 


Data 

censored no. 
total no. 


Cox's prob.(%) 
joint 
z only 


i only 


Redshift, and Optical Luminosity 
EM 

Intercept 
Slope for z 
Slope for i 
Standard Deviation 


Radio QSOs 

P(tot)<0 . 01 

1.1610.03 

11/156 

P(z)<0 . 01 

-0.7310.15 


P(i o )<0.01 

0.2510.03 

0.17 

Optical QSOs 

P(tot)-0 . 01 

1.5610.06 

73/102 

P(z)-53 

0.2110.28 


P(i )-42 
o 

0.0710.06 

0.20 

X-ray AGNs 

P(tot)<0 . 01 

1.2010.02 

0/123 

P(z)<0 . 01 

-1.2910.13 


P(i )<0.01 

0.2710.03 

0.13 

Radio BL Lacs 

P(tot)-0.01 

1 . 1010 . 06 

0/24 

P(z)=0 . 02 

-1.3910.26 


P(i )-0.01 
o 

0.2610.05 

0.15 
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X-ray BL Lacs 
0/10 


P(tot)-23 

P(z)-ll 

P(i 0 )-77 


0 . 90±0 . 12 
-0 . 48±0 . 32 
0 . 08±0 . 90 
0.15 
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Table .9c : Correlation among a , Redshift, and X-ray Luminosity 


Data 

censored no . 
total no . 


ox 

Cox' s prob . (%) 
Joint 
z only 


I only 


EM 

Intercept 
Slope for z 
Slope for $. 


Standard Deviation 


X-ray AGNs 

P(tot)<0 . 01 

4.5310.58 

0/123 

P(z)<0 . 01 

0.8710.20 


P(J )<0.01 

X 

-0.2110.04 



0.15 

BL Lacs 

P(tot)— 51 

3.5211.50 

0/24 

P(z)-76 

0.4710.53 


P(i )-47 

X 

-0.1410.10 



0.21 

X-ray BL Lac 

P(tot)-21 

1.9811.07 

0/10 

P(z)-23 

-0.2210.34 


P(i )-45 

-0.0710.07 


X 

0.10 
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Table 10a : Comparison of Radio Quasars and BL Lac Objects 


Property 


Mean and Standard Deviation 


censored no . 


total no . 


Radio luminosity QSO : <log(i^)> - 33.88±0.128 (46/156) 

BL Lac: <log(i )> - 33.8110.280 (0/24) 

Gehan test : P = 79% 

Logrank test: P - 91% 

Optical luminosity QSO : <log(i Q )> - 31.2610.06 (0/156) 

BL Lac: <log(i o )> = 30.5410.23 (0/24) 

Gehan test : P = 0.8% 

Logrank test: P = 0.01% 

X-ray luminosity QSO : <log(i x )> = 45.6310.06 (11/156) 

BL Lac: <log(Lx)> = 44.8910.24 (0/24) 

Gehan test : P = 0.6% 

Logrank test: P = 0.01% 

a ro QSO : <a ro > ” 0.6410.01 (46/156) 

BL Lac: <* rQ > = 0.6210.02 (2/50) 

Gehan test : P = 20% 

Logrank test: P == 14% 

a QSO : <a > - 1 . 43+0 . 02/-0 . 03 (11/156) 

OX ox 

BL Lac: <a ox > - 1.4610.03 (2/50) 
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Gehan test : P - 29% 


Logrank test: P - 25% 
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Table 10b : Comparison of Radio Quasars and BL Lac Objects (0.8<z<1.7) 


Property 


Redshift 


Radio luminosity 


Optical luminosity 


X-ray luminosity 


a 

ro 


Mean and Standard Deviation 


censored no . 


total no. 


QSO : <z> - 1 . 11±0 . 03 (0/57) 

BL Lac: <z> - 1.08±0.10 (0/8) 

Gehan test : P - 42% 

Logrank test: P - 81% 

QSO : <log(i r )> - 34.40±0. 12 (16/57) 

BL Lac: <log(i r )> - 35.0410.12 (0/8) 


Gehan test : 

P « 3% 


Logrank test: 

P - 0.6% 


QSO : <log(i o )> - 31.3410.06 

(0/57) 

BL Lac: <log(i Q )> - 31.3810.18 

(0/8) 

Gehan test : 

P - 82% 


Logrank test: 

P - 92% 



QSO : <log(i x )> - 45.7810.05 

(3/57) 

BL Lac: <log(^)> - 45.9210.17 

(0/8) 

Gehan test : 

P = 59% 


Logrank test: 

P - 44% 



QSO : <a ro > - 0.6010.02 (16/57) 

BL Lac: <<* ro > - 0.6910.03 (0/8) 
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Gehan test : P - 8% 


Logrank test: P - 3% 

QSO : <a Qx > - 1 . 37±0 . 02 (3/57) 

BL Lac: <a ox > - 1.39+0.08/-0.10 (0/8) 
Gehan test : P - 73% 

Logrank test: P - 60% 
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Figure Captions 


Figure 1(a) Plot of the 0.5-4. 5 keV X-ray luminosity against the 5 GHz radio 
luminosity density of the radio selected quasars. The detected 
points are represented by circles and the upper limits are 
represented by bars in this figure and all following figures. 

The solid line in this and following two figures is the 
regression line log(.2 x )=a+51og(2^) obtained by the application 
of the BHK method. 

1(b) Plot of the 0.5-4. 5 keV X-ray luminosity against the 5 GHz radio 
luminosity density of the optically selected quasars. 

1(c) Plot of the 0.5-4. 5 keV X-ray luminosity against the 5 GHz radio 
luminosity density of the X-ray selected AGNs. 

1(d) Plot of the 0.5-4. 5 keV X-ray luminosity against the 5 GHz radio 
luminosity density of the BL Lac objects. The radio selected BL 
Lac objects are represented by circles and the X-ay selected BL 
Lac objects are represented by boxes. The solid line is the 
regression for the radio selected BL Lac objects and the dashed 
line is the regression line for the X-ray selected BL Lacs. 

Both regression lines are assumed to have a form 
log(i x )-a+&log(^^) and computed by a least square method. 

Figure 2(a) Plot of the 0.5-4. 5 keV X-ray luminosity against 2500 A optical 
luminosity density of the radio selected quasars. The solid line 
in this and the following two figures is the regression line 
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log(i x )-a+blog(i^) by the EM algorithm with normal distribution. 

2(b) Plot of the 0.5-4. 5 keV X-ray luminosity against 2500 A optical 
luminosity density of the optically selected quasars. 

2(c) Plot of the 0.5-4. 5 keV X-ray luminosity against 2500 A optical 
luminosity density of the X-ray selected AGNs. 

2(d) Plot of the 0.5-4. 5 keV X-ray luminosity against 2500 A optical 
luminosity density of the BL Lac objects. The radio selected BL 
Lac objects are represented by circles and the X-ay selected BL 
Lac objects are represented by boxes. The solid line is the 
regression line for the radio selected BL Lac objects and the 
dashed line is the regression line for the X-ray selected BL 
Lacs . 


Figure 3(a) Plot of the 5 GHz radio luminosity density against 2500 A 

optical luminosity density of the radio selected quasars. The 
solid line here and in the following figures is the regression 
line log(f^_)-a+blog(^ o ) by the EM algorithm with the normal 
distribution. 

3(b) Plot of the 5 GHz radio luminosity density against 2500 A 

optical luminosity density of the optically selected quasars. 

The uncertainty of the intercept coefficient is relatively large 
(0.32 ± 0.65). 
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3(c) Plot of the 5 GHz radio luminosity density against 2500 A 


optical luminosity density of the X-ray selected AGNs. The 
uncertainty of the intercept coefficient is relatively large (- 
34.5 ± 12.0) . 

3(d) Plot of the 5 GHz radio luminosity density against 2500 A 

optical luminosity density of the BL Lac objects. The radio 
selected BL Lac objects are represented by circles and the X-ay 
selected BL Lac objects are represented by boxes. The solid line 
is the regression line for the radio selected BL Lac objects and 
the dashed line is the regression line for the X-ray selected BL 
Lacs . 

Figure 4(a) Plot of the against the of the radio selected quasars. 

The solid line is the regression line a =a+bo by the 

ox ro J 

application the BHK method. The bars in this figure and 
following four figures represent lower limits. 

4(b) Plot of the a ^ against the of the optically selected 

quasars. Because of the large uncertainty for the intercept here 
and below, we do not fit a regression line. 

4(c) Plot of the a against the a of the X-ray AGNs. 

ox ° ro J 

4(d) Plot of the against the a of the BL Lac objects. The radio 

selected BL Lac objects are represented by circles and the X-ay 
selected BL Lac objects are represented by boxes. The solid line 
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is the regression line for the radio selected BL Lac objects and 
the dashed line is the regression line for the X-ray selected BL 
Lac objects. 

4(e) Plot of the a ^ against the a rQ of the compact radio sources. 

The arrow in this figure represents an upper limit. 
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